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The asymptotic structure of the Kohn-Sham exchange potential Vx (r) in the classically forbidden 
region of a metal surface is investigated, together with that of the Slater exchange potential Vf (r) 
and those of the approximate Krieger-Li-Iafrate (r) and Harbola-Sahni 144, (r) exchange poten¬ 
tials. Particularly, the former is shown to have the form of Vx(z —^ oo) = —cxxlz with a constant 
dependent only of bulk electron density. The same result in previous work is thus confirmed; in the 
meanwhile controversy raised recently gets resolved. The structure of the exchange hole px{r,r') 
is examined, and the delocalization of it in the metal bulk when the electron is at large distance 
from the metal surface is demonstrated with analytical expressions. The asymptotic structures 
of Vx{r), Vx{r), and M4(r) at a slab metal surface are also investigated. Particularly, 

Vx{z —^ oo) = —Ijz in the slab case. The distinction in this respect between the semi-infinite and 
the slab metal surfaces is elucidated. 

PACS numbers: 71.15.Mb, 71.10.-w, 73.20.-r 


I. INTRODUCTION 


The static (classical) charge potential vanishes expo¬ 
nentially far outside a metal surface. The electronic 
structure in the classically forbidden region as a con¬ 
sequence depends strongly on exchange and correlation 
(xc) effects. Quantum mechanical many-body effects 
thus play a major role in this region. In the seminal work 
of application of the Kohn-Sham (KS) density functional 
theory (pFT) [l[ to the metal surface problem, Lang and 
Kohn [3 assumed that the KS local xc potential Vxciz) 
decays like the classical image potential Vim{z) = — l/4z, 
with z being the distance from the metal surface. (Sim¬ 
ilar assumption had been made earlier in Bardeen’s pi¬ 
oneering study of the electronic properties of the metal 
surface Q.) The KS local xc potential Ua;c(i") is equal 
to the functional derivative of the KS xc energy func¬ 
tional Exc[p] with respect to the electron density p{r): 
Vxc{jt) = 6Exc[p]/Sp{r). It was found that Vxc{z) calcu¬ 
lated by Lang and Kohn in the local density approxima¬ 
tion (LDA) to Exc[p] turned out to have an exponential 
decay at large distance from the metal surface, recognized 
hence as one of shortcomings of the LDA [2. The subject 
of the asymptote of Vxc(z) at the metal surface has there¬ 
after attracted long-standing research interest . It was 
claimed in Ref. and Ref. @ that Vxc{z) ^ ~l/(4-2) 
which was attributed to the correlation potential vdz). 
The exchange potential Vx{z) was claimed to decay ex¬ 
ponentially in the former and as ^ in the latter, 

respectively. 

The exponential and the ~ — I/z^ behavior of the ex¬ 
change potential Vx{z) was questioned in Ref. Q, and 
it was numerically demonstrated that at least the com¬ 
ponent of Vx(z) arising purely from the exchange hole, 
[which is in fact the Harbola-Sahni approximate exchange 
potential Wx(z) discussed later in this paper,] has an 
image-potential-like behavior though possibly not the ex¬ 
act form of Vim{z). The result in Ref. Q was corrobo¬ 
rated by Solomatin and Sahni Q who, based on the inte¬ 


gral equation for the optimized effective potential (OEP) 
0 [also known as the 0PM (optimized potential method) 
in and various literature], analytically showed that 

Vx{z ^ oo) =-ax-, (1) 

z 


with 


/3^ — 1 r ln{l3‘^ — 1) 
2/32 [ “ 7r(/32 - I)i/2j’ 


( 2 ) 


where /3 = sJWjeF, fK the surface-barrier height, and 
tp the bulk Fermi energy. This result had been further 
confirmed in Ref. 0- ' The issue of the asymptote of 
the full KS exchan ge-c orrelation potential Vxc{z) was also 
addressed in Ref. [lfll |. 

Both the calculations in Refs. Q and for Vx{z) at 
large z had been carried out exactly (with no approxima¬ 
tion employed). The one in Ref. [1^ made the use of a 
different method by solving the Dyson equation with the 
use of the exact exchange part of the electron self¬ 
energy. The agreement between Ref. [8| and Ref. [l^ 
on the result shown in Eq. m strongly indicates its cor¬ 
rectness. Recently Horowitz et al. however claimed m 
that asymptotically Vx(z) = — I/z, but later [I^ that 
Vx{z) has an asymptotic form of Vx{z) ^ Inz/z. The 
result in Eq. ([2| was hence challenged. The method 
used in Refs. [ll|) 0 f^e same OEP method used 
previously in Ref. [8|. However there was a subtle but 
distinct difference in the technique between the approach 
in Ref. Q and that in Refs, [ill [13 • The calculation in 
Ref. [3 made the direct use of the integral equation for 
OEP, but those in Refs. [ll|i El ^he use instead of the 
OEP method formulated in terms of differential equa¬ 
tions for the so-called orbital-shifts [cf. Eq. (l2^ below] 
[TMl- There are certain advantages in this formulation 
in that the quantities of the orbital-shifts are compara¬ 
tively amenable to analytical or numerical study. This 
will also get illustrated in this work. Indeed, we attempt 
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to resolve the controversy by carrying out further inves¬ 
tigation based precisely on this formulation of the OEP 
method. To this end, we have managed to establish, 
for limiting large 2 , an identity between Vx{z) and the 
planar-momentum averaged orbital-dependent exchange 
potential Uxk{z) k = kp, where k is the perpendicular 
component of the electron momentum, and kp the bulk 
Fermi momentum. The identity, reported in Eq. (1421) [for 
the case of the slab metal surface see instead that in Eq. 
(ED]. is one of the key results, with the aid of which the 
following development can be made fairly smoothly. The 
quantity Uxkpiz) is then shown to have the asymptotic 
form of Uxkpiz) = —axjz. In this manner, we confirm 
once again the result in Eq. (HD- The controversy raised 
in Ref. [l^ is hence resolved. It is shown in Sec. V 
that one component [termed 14 ®^*-^*(z)] of Vx{z) which 
was ignored in the study in Ref. [l^ actually also makes 
a leading-order contribution to Vx(z). In addition, it is 
shown that Vx(z = 00) = 0 at the semi-infinite metal 
surface, i.e., Eq. (CD is exact to the leading order. Since 
Vxciz = 00) = 0 (which can always be made true), it 
follows that Vc{z = 00) = 0 . 

The Slater exchange potential (r) has been regarded 
also as an approximation to the KS exchange potential. 
In Ref. Q it was shown that Vf(z) [l^ has the following 
asymptotic form at large distance from the metal surface: 

V^{z ^ 00 ) =-2ax^. (3) 

There exists no controversy in the literature over this 
result. We give further verification of it. 

The quite nontrivial asymptotic structure of Vx{z) in 
Eq. (HD was pointed out to arise from the delocaliza¬ 
tion of the exchange hole (also known as the Fermi hole) 
pAr,r')^. It was shown numerically that the exchange 
hole is spread throughout the entire metal bulk when the 
electron is in the classically forbidden region [13 . Ex¬ 
pressions for the exchange hole with the planar-positions 
averaged, valid for arbitrary electron positions, are re¬ 
ported. Especially, the delocalization of the exchange 
hole at the semi-infinite metal surface for limiting large 
electron positions is illustrated with analytical expres¬ 
sions. This might be of help to shed further light on 
the curious exchange effects in the classically forbidden 
region of the metal surface. 

Two of other approximate exchange potentials, the 
Krieger-Li-Iafrate (KLI) V)5^^(r) potentia l fi^ . and 
the Harbola-Sahni (HS) 114 (r) potential [^, are also 
surveyed. They will be introduced in Sec. V. is 

well-known as an extensively employed substitute for the 
OEP. It has the same bulk limit as Vx(z), a fact crucial for 
calculation practice for the metal surface. Vj^^^{z) turns 
out to have an asymptotic form of Vj^^^(z) ^ In zjz and 
hence deviate from Vx{z) in the classically forbidden re¬ 
gion, indicating that an improvement is required there. 
On the other side, the asymptotic behavior of Wx{z) has 
been given a full-fledged study in Ref. [2lj, and it was 
shown that Wx(z) deviates from Vx(z) in both the metal 


bulk and the classically forbidden region, though rela¬ 
tively mildly in the latter in the form of Wx (z) ^ —aw jz, 
where the constant aw is given in Eq. (1551) . 

A great deal of work on the electronic structure at the 
metal surface has been carried out on a jellium metal 
slab instead. The asymptotic behavior of Vx{z) at large 
distance from the slab surface is also examined in this 
work, and it is shown that 

Ua;(2:-t 00) =-i. ( 4 ) 

z 

As mentioned pre viously, the same result was also re¬ 
ported in Ref. [ij. In this sense it gets confirmed here. 
It is necessary to remark that the study in Ref. m 
was performed on the slab metal surface, but clearly also 
aimed at obtaining knowledge equally valuable for the 
semi-infinite metal surface. We, however, point out that 
Eq. (jT]) is valid only for the slab surface and can not 
be naively extrapolated for the semi-infinite surface. In 
fact, the slab surface can virtually be regarded as a finite 
system and as a consequence Eq. (SD can be obtained 
by a direct multipole expansion of the Coulomb inter¬ 
action. Furthermore it is shown that the result in Eq. 
(SD is in fact subject to a possible difference of nonzero 
constant, i.e., Vx{z = 00 ) 4 0 [see Eq. (11031) and the 
discussion below it], though it can be made to vanish in 
an exchange-only calculation via a shift in Vxiz). 

In Sec. II we introduce the OEP method to the metal 
surface problem. The properties of Vx{z), V£{z), and 
the exchange energy density ex{z) are examined in Sec. 
III. The properties of px{r,r') are examined in Sec. IV, 
and those of and II4(z) in Sec. V. The slab 

surface case is considered in Sec. VI. In Appendix A we 
prove a mathematical statement of Eq. (1241) proposed 
in Sec. II. Appendix B contains the calculation for the 
asymptotic form of Uxkpi^)! Appendix C that for V^{z), 
and Appendix D that for the planar-position averages of 
Px{r,r'). The paper is concluded in Sec. VII. 

II. PRELIMINARIES 

In the jellium Q and structureless-pseudopotential 
models of a metal surface with a uniform positive back¬ 
ground of charge 

the KS orbitals are of the form 

( 6 ) 

In Eq. ([S]) (k|| , X||) are the planar components of the 
momentum and position, and (fc, z) the perpendicular 
components, i.e., k = k||-|- kez and r = X||-|- ze^, where 
is the unit vector perpendicular to the metal surface. 
The magnitude of k, on the other hand, will be explicitly 
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denoted with |k| (= + fc^). A and L in Eq. ([S|) de¬ 

note the planar normalization area and the perpendicular 
normalization length, respectively. 


A. Preparatory materials 


The (pkiz) obeys the differential equation 


1 52 
' 2 ^ 


+ w + V(z) Vxc{z) 


(t^k{z) = -k‘^(j)kiz), (7) 


where V(z) is the static (classical) charge potential which 
vanishes exponentially at large distance from the metal 
surface. (/)k{z) has the following asymptotic forms, 


4>kiz) ^ sin[fcz-f d(fc)] for z —>■ — 00 , (8a) 

^ Pk{z)e~'^^ for z —>■ 00 , (8b) 


It follows from the substitution of Eq. ® into Eq. da 
that 


Px{r,r') =- 


1 


dk j dk!e{kF - |k|) 


87r®p(z) 


(15) 


where 


^k,k'{z,z') = (j)l{z)(l)k{z')(i)k'{z)(i)l,{z'). (16) 

Alternatively, the substitution of Eq. da instead yields 


Px{v, r') = - 


nkp pkp 


pKf 

dk / dk'XX'^k,k'{z, z') 

Jo 


t:'^p{z) j 
Ji(A|x|| -xj||)Ji(A'|x|| -xj|| 
Ivn -V^ I 2 


(17) 


where S{k) is the phase shift due to the metal surface, 
K = \/21E - F, and Pk{z) is a power function . 00 
denotes the positive infinity. 

The Dirac density matrix is defined as 7 s(r, r') = 
2^-(/)i(r)(^*(r'). At the metal surface, it has the fol¬ 
lowing form: 


By the use of the identities: 


where A = ^k'f — kJ and X' = \/^f — k'^. 


B. The OEP method 


The exchange-only OEP Q has been proved to be 
equal to the KS exchange potential (r) . We briefly 

outline the method below, starting with an introduction 
to the orbital-dependent exchange potentials [3 [3 


1 ^Ex[{4>j}] 




<^i(r) ,5<^*(r) 


(18) 


I = 2^Jo(g|x|| -x|||); (10) 


dl(A|X|| — Xy I) 

/ *||^||->^o(A:|||x||-x|||) = A— -—, (11) 

Jo 1^11 ^ 11 1 

where Jq and Ji are the zeroth order and the first order 
Bessel functions, respectively, we can rewrite Eq. m as 


2 *^i(^|xii — Xm I) 

ls{r,r') = — dkX(t)k{z)(l)l{z') —^—-( 12 ) 

^ Jo 1^11 ^ 11 1 

The expression for the electron density is 


1 

pi^) = - / 

71-2 Jo 


dkX'^\(j)k{z)\‘^. 


(13) 


[The density is equal to the diagonal component of the 
density matrix, and Eq. m follows directly from Eq. 
m with the aid of the fact that Ji{x) —>■ l/2a: as x —>■ 0.] 
Finally we formulate for the metal surface problem the 
exchange hole which is defined as 


Px{r,r') = 


l7^(rT')P 

2p(r) 


(14) 


and the orbital shifts 


3=Xi 





(19) 


Here (j)i (r) and are the KS orbitals and the correspond¬ 
ing eigenenergies, respectively, and Ex[{4>jY\ i® the ex¬ 
change energy functional of the orbitals. The orbital- 
dependent exchange potentials Uxi{r) in Eq. m are of 
the same form as the Hartree-Fock potentials but con¬ 
structed from the KS orbitals (j)i{r). Explicitly, 


Uxi{r) = 


t f , ns{r,r')(j),{r') 

2Mr)J |r-r'| 


( 20 ) 


The central equation in the OEP method is 


occ. 

i 


where c.c. denotes the complex conjugate of the previous 
term. Equations m and (EU, together with the KS 
equations for the orbitals [cf. Eq. ([T])] [I|, build the 
self-consistent calculation scheme for the exchange-only 
OEP. 

In self-consistent calculation for the OEP, evaluating 
'0i(r) from Eq. (HU) is highly impractical. Equation m 
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therefore is usually rewritten in the form of a differential 
equation [l3l - [l5l| . In our case, 



where 


W + V{z)+ Vxc{z) - -k^ '0k(r) 
[-Va:{z) + u,jk(r) + £>k]d'k(r), 


( 22 ) 


-^k — (^kl'^^a; '^xkl^k)- 


It is demonstrated in Appendix A that Ma;k(r) is inde¬ 
pendent of X||, i.e., it is a function of z only: 


Wxk(r) = Wxk(2)- 


(24) 


The variables X|| and z accordingly can be separated in 
Eq. (|2^ and V'k(r) hence has the form 


V'k(r) = 



(25) 


Correspondingly, with the aid of Eq. ([5]) , it follows that 


+ W+ V{z) + v^c{z)-]^k'^ V'k(2:) 

= [-Vxiz) + Uxkiz) + Dkjipkiz). (26) 


We average over k|| on both sides of the preceding equa¬ 
tion and obtain 


I 

“2^ 


W + V{z) +Vxc{z) 



ipk{z) 


[-Vx{z) + Uxk{z) + Dk]4>k{z). (27) 


Explicitly sm, 


with 


TT \kp 


\ 1 1 - , 11 + a; I 


Solely for later reference, we also list 


Dk 




(33) 


(34) 


(35) 


and note that IIk||k|=feF = 0. Expression for will also 
be used in the later development: 


D 


b 

k 


2kp 

ttA^ , 



1 


(ik||F 



k/kp 


dxxF(x)\ 


(36) 


explicitly: 



—2kp + 


{kp + k)(2kp — k) 

- kfZl ->" 


(kp — k)(2kp -|- /c) ^ kp — k 
kp -|- k 2kp 


kp F k 
2kp 

■ ■ (37) 


C. Identity Vx(z) = Uxkp(z) for z —^ oo 


The planar-momentum averages Uxk (z) of the orbital- 
dependent exchange potentials UxU (r) [rewritten as Uxic 
(z) according to Eq. (El], 



has been introduced into Eq. (|27)) together with the sim¬ 
ilar planar-momentum averages ipk(z) of ipuiz), and Dk 
of i?k. 

We write 

= (<^k|u,|<(.k); 7?k = (<('k|u.k|<(>k). (29) 

Accordingly, it follows from Eq. (1^ that 

Dk = Dl-Dl (30) 

with 

Dt=2{(l>k\vx\M/L-: Dl=2{(l>k\uxk\^k)/L. (31) 

For the semi-infinite metal surface, evidently takes 
the bulk value = —kp/n of Vx(z), (and is accord¬ 

ingly independent of A:.) On the other hand. 


We consider Eq. (071) for limiting large z. On the rhs 
of it the term Uxk(z)4>k(z) which is due to the orbital- 
dependent exchange potential has the asymptotic form of 
fk(z)4>kF(z), with fk(z) a power function. Its existence 
implies that ipkiz) must have the analogous asymptotic 
form: 

'>Pk(z) = gk(z)(j>kF(z), (38) 

where gk(z) is also a power function obeying 

gk(z) - <l}ll(z)<j)'kF{z)g'k(z) - ]^gk(z) 

= + Uxk(z) + Dk\4>k(z). (39) 

The primes denote the derivatives with respect to z. 
Equation El is valid for all k. It turns out that its 
special case a.t k = kp solely is sufficient to determine 
the asymptotic structure of Vx(z). Indeed, 

(/>kF(^)<f>kF(^)3kF(^) + iffkpi^) 

= Vx (z) - Uxkp (z) - Dkp ■ (40) 


D 


b _ 
k — 



(32) 


By using the fact that gkpiz) —t 0, [with the understand¬ 
ing that the homogeneous solution of tpk(z) in Eq. (1771) 


|k-k'|2- 
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is excluded,] one has g'l^iz) ^ o(l/z) for large z. 

On the other hand, it will be shown later in Eq. (|T7l) that 
Uxkpiz) ^ 0(l/z). One hence obtains, from Eq. (l40l) . 
for limiting large z 

Vx{z) = Uxkpiz) + Dkj,. (41) 

Since = —kp/ir, and = —kp/TT according to Eq. 
dSzl), one has correspondingly, from Eq. dSQl), Dk^ = 0 
for the semi-infinite metal surface. This leads to one of 
the central identities: 


Vx{z) = Uxkpiz), for z^oo. (42) 


We remarked that D^p 7 ^ 0 for the slab metal surface, 
and we then have Eq. (HD) only. 

Equation (1421) suggests that knowledge about the 
quantity Uxkp (z) could be useful in carrying through the 
reminder of our task. In this regard we include here some 
expressions for it. In general, for k ^ kp and any well- 
behaved function /(k) = /(k||, k) [uxiaiz) in the present 
case], 


/^dk||/(k) = 7rAV(k|| =0,A: = fc^). 

^0 


Hence from Eq. ()28)) . it follows that 

'^xkp (-2^) |k|| =^0,fc—fcjpr • 

Accordingly from Eqs. dsni) and (l24l) one has 


^xkp{z^ — 


1 


‘^4^kp (■^) 


dr' 


,7s(r,r')(()fep(z') 


r - r' 


(43) 

(44) 

(45) 


III. ASYMPTOTIC STRUCTURES OF Vx{z) 
AND Uf(z) 


which, together with Eq. (|42|) . leads to Eq. ([T]), fulfilling 
the main object of this work. (The controversy raised in 
Ref. [H will get further elucidated in Sec. V.) Equation 
(H71) which is exact at limiting large z, together with Eq. 
( 1 ^ . leads also to the conclusion that Ua;(oo) = 0 at the 
semi-infinite metal surface. In other words, Eq. CD is 
exact to the leading order. We assume that Vxc{oo) = 0 
at the metal surface. Accordingly we also have Uc(oo) 
= 0. Eurther insight into these facts will be provided in 
Sec. VI in connection with those for the slab case. 

The Slater exchange potential V£ (r) [l^ is defined as 

r/(r) = y drVx(r,r')|^^^. (48) 

For our metal surface problem we substitute Eq. (1171) for 
Px (r,r0 iiito Eq. (l4^ . It follows that 


Vx^(z) = 


A nkp pkp 

- 3 / X / dk(f)k{z) / dk'(j3l,{z)XX' 

P\^) Jo Jo 

POO 

dz'(j)l{z')(j)k'{z') 


dr'u 


0 r'^(z-z02 + r '2 


Ji(Arf|)Ji(AV[|). 

(49) 


The asymptotic structure of Vjf (z) at large distance from 
the metal surface is calculated in Appendix C and the 
final result is shown in Eq. m i- Correspondingly, 
the exchange energy density per unit volume ea;(r) = 
4p(r)V)f (r) has the asymptotic form 0: 

ea;(z-7 00 ) =(50) 

z 

We note that 


Quantities similar to Uxkp (z) have been calculated in 
Refs. We apply a different approach to per¬ 

form the calculation, which can serve also as a kind of 
verification. To this end, we substitute Eq. (HU) into Eq. 
(HSl) . It follows that 


Uxkp (^) — 


fkp 


T^ 4 >kp {z) Jq 

dz' (t)l{z')(t)kp{z' 


dkX(j)k{z) 


dr'n 


0 ^(z - z'Y + rj 


= Ji(Ar:,). (46) 


Alternatively Eq. (l46ll can also be obtained from the 
substitution of Eq. (IA3I) into Eq. (H^ . The aid of the 
fact that Jo(0) = I is needed then. 

The asymptotic form of Uxkp{z) is examined in Ap¬ 
pendix B with the following result. 


Uxkp{z -)■ 00 ) 


I 


— 


(47) 


p{z^oo) = ^§^j^\(t)kp{z)\'^. (51) 

It is well known that in a finite system r 7 (r)has a 
long-range form of Vx{r —)> 00 ) = — I/r [ 1 , 0 , IH, 1 ^ . 
(For simplicity, only a spherically symmetric system is 
discussed here. For recent progress made on this sub¬ 
ject, see Refs. Bin and the discussions later in Sec. 
VI.) The Slater exchange potential has exactly the same 
long-range form of Vjf(r —7 00 ) = —I/r. The identical 
long-range form of Vx (r) and U/ (r) results from the fact 
that in the finite system the exchange hole is well local¬ 
ized near the system. The asymptotic structure of Vxir) 
is essentially determined by that of the orbital-dependent 
exchange potential U xm(r) o f the highest occupied orbital 
(denoted by m) iS [ 2 ^ [ 2 ^ . The asymptotic structure of 
both Mxm(r) and Vjf(r) in turn can be attributed to the 
(orbital-dependent) exchange hole in terms of the lead¬ 
ing order contribution of the multipole expansion of the 
Coulomb interaction. For the semi-infinite metal surface 
this is no more true since the exchange hole is delocal¬ 
ized and spread throughout the entire bulk region [13, 
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and as a consequence the multipole expansion argument 
does not apply anymore. This explains the difference in 
the asymptotic structures of Vx(z) and V^{z) at large 
z from the metal surface, as shown in Eq. © and Eq. 

respectively. The delocalization of the exchange hole 
at large electron positions will be examined in the next 
section. 


IV. DELOCALIZATION OF px(r,r') 


We take planar-position average of px{r,r'): 


Pxiz,z') 



dx||P3,(r,r'). 


(52) 


The information about px{z,z') clearly is sufficient to 
serve the present purpose. By the substitution of Eq. 
(fTSl) into Eq. (15^ . Px{z, z') can be shown as 


We substitute Eq. (ITSl) into Eq. (l56t and write bx{z,z') 
as 


1 


Sir^piz) Jo 

J dkj dk'eikp - \k\)e{kF - |k'|) 

^k,k'{z, ^ 


(57) 


By introducing the transforms: q = k|| — kj| and K = 
(k||-fkj| )/2, one may rewrite bx{z, z') as 


bx{z,z') =- 


1 


r^F 


/ kp 

dk'^k,k'{z, z') 

-kp 


Sn^piz) J_kp J-kp 

poo 


J dqFiq) d|x||-X|||e 
where the function F{q) is defined as 


Px(z,z')= - 


nkp 


dky 


TT^piz) Jo 

'k 

dk'[(kk,k'{z,z') + $fe,fe/(z',z)]. (53) 


Equation (15^ is valid for all electroii positions z. This 
expression was first reported in Ref. [T^ and the reader 
is referred to that work for detailed derivations. It is not 
difficult to see that the two terms in Eq. (155)) make the 
same contribution, and hence 


Px{z,z') = - 


f^kp 


Tr^p{z) 


dkX^ (j)l{z)(l)k{z') 


dk'(j)k'{z)(j)*k,{z'). (54) 


F{q) = y'dK0(A-|K-f q/2|)0(A'-|K-q/2|). (59) 

By using Eq. (flUl) and further the following identity: 

f°° 1 

/ c?|x||-xj||Jo( 9 |x||-xj||) =-, (60) 

Jo d 


Eq. 


may be further rewritten as 


bx{z,z')= - 


I 


pkp 


7r®p(z) 

‘kp 


dk 


poo 

dk'‘Pk,k'{z,z') / dqF{q). (61) 
Jo 


Derivation based on Eq. (I5T]) . detailed in Appendix D, 
yields the following result for bx{z, z'): 


We are interested especially in the asymptotic struc¬ 
ture of px{z,z'). It is demonstrated in Appendix D that 

Px{z^oo,z') = — --IJikpiz')]'^. (55) 

TTC Z 

Equation (1551) displays remarkably the delocalization of 
the exchange hole in the metal bulk. [The reader is re¬ 
ferred to Eq. (l8al) for the behavior of Jikp {z) in the deep 
bulk region.] In the meanwhile, the amplitude of Pxiz, z') 
decays as ^ I/z with the electron position z. 

Both planar positions xy and xj| of the electron and the 
hole have been averaged, respectively, in Eq. (1521) . But 
evidently px{r,r') is a function of xy—x|| only, rather 
than of X|| and xj| separately. Another type of average 
over the planar positions therefore might be equally ca¬ 
pable of featuring the exchange hole at the metal surface 
which is defined in the following manner (T7j |. 


bxiz,z') 


16 

3'!r^p{z) 


[G(z, z') -I- G(z', 


^)], 


(62) 


where 

G{z,z') 


pkp pkp 

/ dk dk'^k,k'{z,z') 
lo Jk 


K 


(A'^_A^) + E(y )(A^ + A'^) 


(63) 


K and E are the complete elliptic integrals of the first 
and the second kinds, respectively [^ : 


K{t) = 



(I — — cc^) 


dx, 


(64) 


poo 

bx;{z,z')= d|x|| - xj||p,^(r,r'). 

^0 


(56) 




I — t^x'^ 


dx. 


(65) 
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We note that so far all results for bx{z, z') hold in general 
for arbitrary position z. Particularly, Eq. (EH) evolves 
into the following asymptotic form for limiting large elec¬ 
tron position z, (see again Appendix D for the demon¬ 
stration) 

b^{z 00, z') = (66) 

where 

x[K(Vy)(2/-l) + E(v/^)(2/ + l)]. (67) 

Equation (1551) displays similar delocalization nature of 
the exchange hole. The amplitude of bx{z,z') decays as 
^ Xjz^l"^^ a little faster than that of px(z^z') shown in 

Eq. EH- 


20^2 

7^5/2 


V. ASYMPTOTIC STRUCTURE OF 
APPROXIMATE EXCHANGE POTENTIALS 


The Slater potential V)f (r) was proposed as an approx¬ 
imation to the orbital-dependent Hartree-Fock exchange 
potentials [l^. Upon the arrival of the KS-DFT, V)f(r) 
was regarded naturally also as an approximation to the 
KS exchange potential Ua;(r). The large-distance struc¬ 
ture of it at the metal surface is shown in Eq. © and 
discussed in Sec. III. In this section we survey those of 
two of other approximate exchange potentials, the KLI 
potential E/"^|(r) and the HS potential 114, (r). The for¬ 
mer reads [ll, fl^ 

V^^\v) = U/(r) + U,^(r), 

where 


( 68 ) 


W = ^ 1] \(j)i{r)\‘^{(l>i\vx - Uxi\(l)i). 

' i 

In passing we mention that [13 

vx{r) = U/^^(r) + U;"*/‘(r), 
and the component I4f^*^*(r), defined as 


(69) 


(70) 


occ. 


r\J ^ 

-V'.(r)v'</':(!■)], (71) 

is sacrificed for the purpose of less calculation labor. [For 
convenience, the same symbols V^{v) and (r) in 

Ref. [l2| are adopted here.] 

We first make a digression to comment on some bulk 
properties of ipk{z) and Vj^^\z). The expression for 
at the metal surface is 


V' 
’ 1 


shift 


(z) = 


f'kp 


2Tr^p{z) Jo 


dky 


J^k{z)^J’k{z) - 4fe(z)^^fc(z) 


.(72) 


On the other hand, in view of Eq. EH, one has 


Vx[z = -oo) = Dl 



(73) 


and 


Uxk{z = -oo) = D\. (74) 

It follows that the rhs of Eq. EH vanishes. With the 
understanding that the homogeneous solution of ipkiz) 
in Eq. is excluded, one thus has 

tpk{z) = 0 for z = —oo. (75) 

Therefore Vf^'^^*'{z), according to Eq. dzH, vanishes in 
the metal bulk: 


= 0 for z = -oo, (76) 

and consequently according to Eq. dlQl), has 

the merit of possessing the same bulk limit as Vx{z): 

= for z = -oo. (77) 

TT 

We next return to the issue of the asymptotic behavior 
of Vj^^^{z) in the classically forbidden region. Vj^^^{z) 
is shown below to deviate strongly from Vx{z) at large 
distance from the metal surface. To this end, we first 
write the expression for Vf{z): 

Vx^(z) = r dk\^\d,k{z)\^Dk. (78) 

^ P\^) Jo 

For limiting large z, Vf{z) turns out to have the form: 


Vfiz^oo) 


27r z 


In 


/ kpz \ 

V^/^y 


-bC'-b21n2- 1 


(79) 


with C = 0.577215 the Euler constant. In obtaining Eq. 
dZH, the following fact 


Dk = —{kp - k) 
47r 


1 - 21n 


kp — k 


2kp 


+ o{kp — k), (80) 


for small kp — k, [which follows from Eq. (1551) ]. has been 
employed. Equation (1751) had also been reported in Ref. 
[l2j ]. Finally, Eq. (155|) - together with Eqs. dH and (1751) . 
yields 


U, 


KLI 


(z —>■ oo) = 


— 1\. f kpz 


[h 


27r z^ V - 1 
-bC'-b21n2-l -2ax 


(81) 


which has the leading form of 0(lnz/z), and hence devi¬ 
ates from Vx{z). 























The HS potential Wx(r) is defined as the work 
done in the Pauli field £x{y) [W^i^) hence also known 
as the Pauli potential in the quantal density functional 
theory (Q-DFT)] [1^, 


Wx{y) = - [ £x{r')-dl', 
J oo 


and 


£x{'r) 



I-') V 



(82) 


(83) 


to Vx{z) [l^. On the contrary, we find that in fact 


14 "'**/* (z 


^)= - 




In 


27r 

+C' + 21n2- 1 


kpz 


1 


+ Otx 


( 88 ) 


The leading order of 0(lnz/z) therefore exactly cancels 
out in Vx{z), resulting in a faster decay of 0(l/z). The 
discrepancy of from Vx{z) in the classically for¬ 

bidden region is also due to the ignored contribution [cf. 
Eq. dZOl)] from V;"'**/*(z). 


The asymptotic structure of W4(z) has been extensively 
investigated in Ref. [2l|, and one of the main results is 

Wx{z) = -awK (84) 

where 

^2_ir^2_2 2 

I 

It might be helpful to mention the relation between 
144(r) and Vxir) [13]: 

u,,(r) = 414 (r) - VE/c^(r), (86) 

where Wf^\r) is the lowest-order correlation-kinetic 
component of Ua;(r). 414 (r) bas been well recognized as 
the component of Ux(r) which arises purely from the ex¬ 
change hole px{r,r') [23,[33|. Clearly 114*-^^(z) must also 
be long-ranged at the metal surface. This is in sharp 
contrast to its behavior in finite systems such as atoms 
and molecules in which (r) decays in a rather short- 
ranged form of ^ 1 /r® [Ig . We wish to further mention 
that it was shown in Ref. [2l| that 414 (z) deviates from 
Vx{z) also in the metal bulk due to the contribution from 
44tc^(z)- Only in the LDA, 144(z) has the same bulk 
value of —kp/iT as Vx{z), as shown in Ref. [3l| |. In this 
connection we list the long-known fact for (z): 



Vf{z} = 


3kp 

27r 


for z = —oo. 


(87) 


VI. ASYMPTOTIC STRUCTURES OF Vx{z), 
vi{z), Uf"(2)> and Wx{z) at the slab 
SURFACE 

So far we have considered only the metal surface with 
the semi-infinite geometry. A great deal of work of the 
electronic structure at the metal surface has been carried 
out on a jellium slab instead. We shall consider the slab 
case in this section. To this end, we first mention two 
studies in which the asymptote of the xc potential in the 
classically forbidden region was addressed as one of the 
key issues. One was reported in Ref. m , mentioned 
previously in the Introduction, in which it was claimed 
that Vx{z —>■ oo) = — I/z asymptotically. Numerical cal¬ 
culation in Ref. [s^ based on the GW approximation to 
the electron exchange-correlation self-energy T^xc, how¬ 
ever, yielded the different result that Vx{z) ^ —Xjz^ and 
Vc{z) ^ —I/(4z). Both studies were performed on the 
slab surface. Part of our effort in this section will be de¬ 
voted to sheding some light on these results. In fact it 
is found that the asymptotic behavior of Vx{z), V£{z), 
VKLi^z)^ and 414 (z) all depends critically on the width 
of the slab. 

We consider a metal slab with a typical width not ex¬ 
ceedingly larger than \p or, in other words, comparable 
to or smaller than \p where Aj’ is the bulk Fermi wave¬ 
length. In this case the discreteness of the eigenenergies 
of the electron in the slab must be taken into account. 
The system can be regarded virtually as a finite one, 
and the well-known conclusion for the finite system that 
Vx{y) ^ —1/r mentioned in Sec. Ill therefore holds. As 
a matter of fact, this already explains the result in Eq. 
O- More explicitly, the Dirac density matrix has the 
following well-known asymptotic form: 


Therefore, strictly speaking, only, among all the 

approximate exchange potentials considered here, has the 
remarkable property with the same bulk value as that of 

Vx{z). 

Finally we would make several comments on the results 
reported in Ref. [ll|. The large-distance structure of 
Vx{z) was attributed solely to Vj^{z) in Ref. [l3|, and 
consequently the conclusion Vx{z) ^ Inz/z was reached. 
On the other side, was claimed to decay as 

^ Inz/z^, and hence make no leading order contribution 


74r,r') = 2(/)™(r)(()4(r') for z^oo. (89) 

From Eq. (unD, by the use of multipole expansion ar¬ 
gument for the Coulomb interaction which evidently is 
applicable here, one immediately has 

Uxra{z ^ oo) = . (90) 

z 

The notation m clearly has the same meaning as kp but is 
used instead to emphasize the discreteness of the eigenen¬ 
ergies. Equation (j?]) then follows from Eq. (|^^ . 
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Additional care is needed in the above discussion for 
Vx{z) in that it is not fully rigorous since actually D^p ^ 
0 in the case of the slab metal surface. Equation (|42ll 
therefore is not valid and we in effect must resort to Eq. 
(1411) instead. This point will get further refined near the 
end of this section [cf. Eq. (|103l) below]. 

In any way, the fact that the planar freedom re¬ 
mains extending to infinity is essentially irrelevant in the 
present discussion. The key point is the dominant be¬ 
havior of the highest occupied orbital and localization 
of the exchange hole near the finite system. The slab 
clearly fails to catch the continuity feature of the elec¬ 
tron eigenenergies which is essential to the semi-infinite 
metal surface. Equation (|4]) thus has limitation in that it 
is valid only for the slab surface and can not be naively 
extrapolated for the semi-infinite surface. Nevertheless, 
a numerical demonstration of Eq. ([4]) for the slab surface 
like that in Ref. m is still quite valuable. 

The effects of the delocalization of the exchange hole 
are negligible in the case of the slab surface. V^{z) con¬ 
sequently has the same asymptotic form as Vx(z): 

E/(z^oo) = -i (91) 

z 

at limiting large z. Correspondingly e^, (z —>■ oo) = 
—p{z)l2z. On the other hand p{z) has the following 
asymptotic form: 

p{z ^ <X)) = 2\(t>kp{z)\'^. (92) 

The reader is referred to Eq. dSB for comparison. Equa¬ 
tion (IM|) can be readily obtained from Eq. (l4^ via keep¬ 
ing the leading-order term in the multipole expansion for 
the Coulomb interaction. It is basically just another il¬ 
lustration of the well-known result V^{r —>• oo) = — 1/r 
for a finite system. Indeed, it follows from Eqs. (11^1) and 
(15^ that, for z —>■ 00 , 


In summary Eq. dsni) [or Eq. (pTI) ] can be obtained for 
the slab case via the argument for the multipole expan¬ 
sion of the Coulomb interaction. The metal-surface fea¬ 
ture of the slab plays no crucial role at this point. There¬ 
fore it is not necessary to resort to detailed derivations 
based on Eq. (ITO [or Eq. . [Equations (HS1) and 
hold for both the cases of the semi-infinite and the slab 
metal surfaces.] Nevertheless such derivations turn out 
to be amazingly simple (due of course also to the feature 
of the finiteness of the system). In the meanwhile they 
could be fairly illuminating. We therefore include one in 
the following mainly for the purpose of illustration. To 
this end, we first copy Eq. (|46l) below with appropriate 
modifications in the form for the slab case: 


'^xkp ('^) — 


Lcjikpiz) 


'^Mkiz) 


dz'<j)l{z')<j)kp{z') 


f-L 

poo 


dr'n 


1 


(z -z'y+ r'2 


Ji(Ar'). (97) 


Only the main domain {—L < z' < 0) has been taken into 
account for the integral over z', which is clearly justified. 
Since z >> z', the denominator ^(z — z'Y + can be 

readily replaced by y^z2”+r^. The dominant contribu¬ 
tion to the integration over k arises from the region of 
kp — k ^ 1/z, and accordingly A >> kp — k. On the 
other hand. The dominant contribution to the integra¬ 
tion over r|| arises from the region of 0 < rJi < 0(A“^) 
since Ji(oo) = 0. Correspondingly ry << z and hence 

^z"^ + r'i^ can be further replaced by z. The integration 
over r|| then turns out simply to be 


Px(r,r') = |(^^(r')|2. (93) 

Equation (PH) immediately follows from Eqs. (pS)) and 
(1^ . Similarly, Eq. (1551) together with Eq. (IM|) yields 


dr[| Ji(Ar||) 


1 

A' 


Equation (l97l) as a consequence becomes 


(98) 


£x{r) = —(94) 

and it then follows from Eq. (l82)) that 

W,{z) = (95) 

z 

possesses the same asymptotic form [but see 
also Eq. (IIOII) ]. Only in this case 14®^*^‘(z), which now 
decays exponentially at large z, is much smaller than 
v^{z) which in contrast approaches a nonzero constant: 

= dJrn for z —>• oo. (96) 

We note once again that kp has the same meaning as m. 
The fact that 7Ik||k|=fc_F = Dkp guarantees no ambiguity 
in the meaning of Dm- 


Mxfcir (-2 -?■ oo) = - 


zL4>kp{z) 




dz'(l>l{z')(l)kp{z'). (99) 


J-L 

One then employs the following equation: 

pO 


J dz'(j)l{z')(j)kp{z') = ^5k,kp- 


( 100 ) 


Equation (15(11) follows (with m equivalent to kp). Equa¬ 
tion (PH) can be obtained from Eq. (P5|) in a similar 
manner. 

Notice that, strictly speaking. 


V^^\z = Dm 


1 


( 101 ) 
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where the term of the constant Dm arises from the (z) 
component on the rhs of Eq. (IS51) . as shown in Eq. (IMl) . 
Furthermore, in accordance to the claims made for finite 
systems in Ref. [l^, such type of constant might also 
possibly occur in Vx{z), since as mentioned above the slab 
can be regarded essentially as a finite system. Indeed, 
since Dkp ^ 0 for the slab surface, one should resort to 
Eq. (I41II instead of Eq. (1421) . Equation m in fact can 
be understood as one of the special cases of the following 
general result: 

Vxir) = Uxmir) + Dm, ( 102 ) 

proposed for the finite system in Ref. [l^. In the light 
ofEq. (I90l) we consequently have 

Va:{z ^ OO) = Dm - l/z, (103) 

instead of Eq. Therefore, strictly speaking, Va;(oo) 
^ 0 now, [see also the discussions in the next para¬ 
graph.] Equation (II03I) can be alternatively obtained 
from Eq. dZOD and Eq. (IIOII) . together with the fact 
that —)> oo) vanishes exponentially. 

Finally we remark that V[j^(oo) = 0 in any case. 
On the other hand, Eq. (1821) automatically gurantees 
Wx(po) = 0. Furthermore, Dm can be made to vanish 
so that Vx{oo) = 0 by shifting ^^(r) in the finite sys¬ 
tem [T^ or the slab of the present case. It is exactly in 
this sense that we justify the result in Eq. ®. However 
we wish also to mention that it remains unclear whether 
Vx{oo) can always be simultaneously shifted to be zero 
in the exact (full) KS scheme in which one has already 
adopted Vxc{oo) = 0. The remarkable fact is that Dm in 
Eq. (]I03p evolves into Hk at jkj = kp of Eq. ((551) and 
vanishes for the semi-infinite metal surface. Thus for the 
semi-infinite metal surface one has definitely v^cioo) = 0 
and Vx{oo) = 0. 

VII. CONCLUSIONS 

By the use of the OEP method, we have established an 
identity between the planar-momentum averaged orbital- 
dependent exchange potential Uxkp{z) and the KS ex¬ 
change potential Vx(z) in the classically forbidden region 
of the metal surface. Based on it, the asymptotic form of 
Vx(z) at large distance from the metal surface has been 
investigated. The result is Vx{z —>■ oo) = —ax/z, which 
resolves the controversy raised recently in the literature. 
The point that Vxc{oo) = 0 and Vx(oo) = 0 hold si¬ 
multaneously gets emphasized and carefully elucidated. 
The asymptotic form of the Slater exchange potential 
—>■ oo) = —2axlz is also verified. The result for 
Vx{z) and that for U/(z) were initially proposed in Ref. 
Q , and the former was verified subsequently in Ref. M- 
The further confirmation in the present work indicates 
beyond doubt that they are correct and the issue is fi¬ 
nally settled. Furthermore, the structure of the exchange 
hole has been examined, and especially the delocalization 


nature of it for an electron far outside the metal surface 
has been demonstrated. It is exactly such delocalization 
that gives rise to the quite nontrivial asymptotic behav¬ 
ior of Vx{z). In addition, the asymptotic structure of the 
approximate KLI exchange potential V^^^{z) and HS 
exchange potential Wx{z) at large z has also been sur¬ 
veyed, which are of the forms Vj^^\z —>■ oo) ^ Inz/z 
and Wx(z oo) ^ —aw jz, respectively. 


As mentioned in the Introduction, common wisdom 
favors the belief that the full Kohn-Sham exchange- 
correlation potential Vxc{z) decays like the classical image 
potential at large distance from the metal surface. Doubt 
however has been casted on it in Ref. [a in which the 
asym ptot e of Vc(z) was also studied. It was shown in 
Ref. [i3 that the asymptotic form of Vxc{z) is not the 
same as the classical image potential. Recent progress 
on the asymptote of Vc(z) has also been made in Ref. 
[ 3 ^ . Unfortunately approximations have been employed 
in the calculations for Vc{z) in both of Refs. a and [3^ 
and it is not clear whether they are fully justified. The 
subject of the asymptote of Vxciz) thus remains not fully 
settled. 


The asymptotic behavior of the exchange potential at 
the metal slab surface has also been investigated. It 
is shown that asymptotically Vx{z), as well as V^{z), 
yKLi^z)^ and Wx{z), depends critically on the width of 
the slab. In particular, if the width is comparable to or 
smaller than Af, the slab can be essentially regarded as 
a finite system and asymptotically Vx{z) = —1/z. Vjf (z), 
y-KLi^z)^ and Wx{z) all have this same form. The ex¬ 
change energy density ex{z) correspondingly approaches 
asymptotically —p{z)j2z. All these facts are inherently 
due to the localization of the exchange hole in the fi¬ 
nite system. While by definition Uaf(z) and Wx(z) van¬ 
ish in the classically forbidden region of both the semi¬ 
infinite and the slab metal surfaces, a careful analysis 
reveals that the —1/z term commences only to the sec¬ 
ond order contribution to Vx(z) and Vj5^^(z) and actu¬ 
ally both Vx(z —>■ 00 ) and V^^^{z —>■ 00 ) approach a 
nonzero constant, i.e., Vx{oo) = V^^\oo) = Dm ^ 0. 
The constant Dm however can always be made to vanish 
in th exchange-only self-consistent calculations by shift¬ 
ing Vx(z). It is precisely in this sense that one has the 
result of Eq. ®. 


This work was supported by the National Science 
Foundation of China. 
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Appendix A: Proof for Eq. (1241) 

We substitute Eq. (fT^ . together with Eq. ([5]), into 
Eq. (l20ll . This results in 


Uxk{r) = - 


7r20k(r) V AL Jq 

/ OO 

dz' (j)l,{z')(j)k{z') 

-OO 


dk'X'(l)k>iz) 


, Ji(A'|x|| - Xm| 

. 


(Al) 


Making the transform rj| = X|| — x|| in the above equation 
and making the use of Eq. m once again, we establish 
Eq. (1^^ explicitly with the following expression: 


Uxkiz) = - 


f^kp 


^4>kiz 


dk'X'(j)k' (z) 


dz'<pt,(z')4>k(z') 

Ji(A'r') 




r'J{z-z'r + 


.'2 


(A2) 


Further algebra then yields 

2 


Uxk(z) = - 


dk'X'4>k' {z) 


T^(t)k{z) Jo 

/ OO 

dz'(j)l,{z)(j)k{z) 

-OO 

, Jo(fc||r[|)Ji(A'rf|) 


> —OO 
^OO 


drn 


(A3) 


l{z-z'Y+r'l 

Appendix B: Derivation for Eq. (1471) 

It is not difficult to see that, for limiting large z, the 
leading contribution to the integral over k on the rhs of 
Eq. (l46l) arises from the region kp — k^ Ijz, and that to 
the integral over z', on the other side, from the region of 
the metal bulk. Accordingly, A >> {kp — k) and hence 
|z — z'l >> 1/A. Since Ji(oo) = 0 and accordingly the 
dominant contribution to the integral over ry arises from 
the region of 0 < ry < (9(A“^), one has |z — z'\ » ry in 
the integral over ry in Eq. (H51) . Thus it follows that 


dr'n 


^(z-zO^+rf 


Ji(Ar[|) 


1^-A| Jo 

1 

A|z-zr 


dr\\Ji{Xr\^) 


Correspondingly, 

Uxkpiz: ^ oo) =- 


rkp 


dk(j>k{z) 


'^J’kp (z) Jo 

1 

/ dz'———(i>l{z')4)kp[z).{m) 

J— OO IZ' Z I 


We next define 


r°° 1 

0fc,fc/(z) = 2/ dz' _ (/)fc,(zO/>fc(^:')- (B3) 

J—OO \Z Z I 

Since, as just mentioned, the integral over z' arises 
mainly from the bulk region, one can make the use of 
Eq. (I8a1) for the orbitals in the above expression and it 
follows that 


dz' 


^k,k'{z)= / , [cos(/c_z'+ (5_) 

J_oo 2^-2:' 

— cos(fc+z'+ (5+)], (B4) 

where k± = k± k', 6± = 6{k) ± 6{k'), and —d stands for 
a negative position near the surface whose exact value is 
irrelevant for z >> d. Indeed, 

Qk,k'{z)= / — [cos{fc_(z - z') + (5_)} 

J —OO ^ 

- cos{(fc+(z - z') + <5+)}],(B5) 

and the upper limit of the integral: z + d, can be readily 
replaced by z. Thus one has 

/■°° dz' 

Qkk’{z) =cos{k-z + 6-) / —-cos(fc_z') 

Jz Z' 

r°° dz' 

+ sin(fc_z + i5_) / —-sin(fc_z') 

Jz z' 

/■“ dz' 

— cos(fc+z + (5+) / —-cos(A:+z') 

Jz z' 

dz' 

— sin(fc-|_z + d-i-) / — sin(fc+z'), (B6) 

J Z ^ 

or 

Qk,k'{z)= — COs(fc_Z + (5_)cz(fc_z) 

— sin(A:_z + S-)si{k-z) 

+ cos(fc+z + 6+)ci{k+z) 

+ sm{k+z + 6+)si{kpz), (B7) 

where si and ci are the sine integral and cosine integral, 
respectively (^ : 



si{x) = — 

r dt . 

/ —smt, 

Jx t 

(B8) 

(Bl) 

ci{x) = — 

/■°° dt 

/ — cost. 

Jx t 

(B9) 
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For z —> 00 the third and fourth terms on the rhs of Eq. 
(1571) will be demonstrated at the end of this appendix 
to make only higher-order contribution to Uxkpi.^)- Con¬ 
sequently these terms can be ignored. Hence one comes 
to 


Qk,k'iz^oo)= — cos{k-z + S-)ci(k-z) 

— sm{k-z + 6-)si{k-z), (BIO) 

We now substitute Eqs. (IB3I) and (IBIOI) (with k replaced 
by kp, and k' by k) into Eq. (IB2() . and obtain 


Uxkpiz^oo) = 




dk(j>k{z) 


T^<Pkp (^) Jo 
cos(a + S-)ci(a) + sin(a -|- d_)si(a)], 

(Bll) 


where S- = S{kp) — S{k), and a = (kp — k)z. Again, as 
mentioned above, for large z the dominant contribution 
to the integral in the above equation arises from the re¬ 
gion kp — k ^ Ijz. Correspondinly —>■ 0 and 4)k{z) 
= where c = — 1. Therefore we fi¬ 

nally obtain Eq. (gZl), with 


1 /■“ 

ax = -/ dae “[ci(a) cos a-I-si(a) sin a]. (B12) 

^ Jo 

The value of ax given in Eq. @ may be obtained via 
carrying out the integral in Eq. (IB12I) (^ : 



dae [ci{a) cos a -I- si (a) sin a] = 


TT -|- 2c In c 

2(1-fc2) ■ 


(B13) 


Solely to draw a connection with the calculations per¬ 
formed in an alternate approach in Ref. , we mention 
the following identity 



du^ 


u 

+ 


—ci{a) cos a — si{a) sin a. (B14) 


Were the third and fourth terms in Eq. (IFH) taken 
into account, they would have made the following contri¬ 
bution to Uxkp (z)' 


-[ cos{2fcirz -I- 6{kp)}ci{2kpz) 

cnz 

+ sin{2kpz + S{kp)}si{2kpz)]. (B15) 


Since ci{2kpz) —?> 0 and si{2kpz) —>■ 0 as z —>■ oo, this 
contribution is therefore of the order o(l/z) and the ne¬ 
glect of these terms in Eq. (lB7l) is thus justified. 


Appendix C: Verification of Eq. ([3]) 

The 4>kiz) is in fact real, and hence the integrals over 
k and k' in Eq. (PI) are clearly symmetric. Accordingly 


we change the domain of the k' integration in Eq. (|4^ 
to fg dk' combined with a corresponding double. Similar 
arguments to those at the beginning of Appendix B then 
lead, for z —>■ 00 , to 



dr'u 


r\sJ{z-z'Y+r'‘ 




1 1 

- -- / dt-Ji(t)Ji(AT/A) 

\z ^ I Jo ^ 

A 

2A'|z-z'r 


(Cl) 


Thus it follows from Eq. (l4^ that, for z —>■ oo. 


= - 


4 ^ 

/ dk(j)k{z)X^ / dk'(l)l,{z) 
P\^> Jo Jo 

nOO 

/ dz'-^4>W)<i^k{z'). (C2) 

J —oo ^ I 


Making the use of Eq. (IB3p in the preceding equation 
yields, for z —?> oo, further 

Vx{z) = -/ dk(j)k{z)\^ 

T^'^p{z) Jo 

dk'c^l,[z)Qk',k{^). (C3) 

Clearly the third and fourth terms in Eq. (jB7p for 
&k>,kiz) can be ignored once again. Accordingly Eq. 
(IBIOI) for <dk',k{z) is substituted into Eq. (IC3I) instead. 
We make further the transform of the integral variable 
k' = k — bjz. It follows then that, for z —>■ oo. 



v!{z) = 




2 \2 


Tr^p{z) 


dk\4)kiz)\ X 


dbe ™[ci(b)cosb + si{b)smb]. (C4) 


We next compare the above equation with Eq. (IB12I) to 
obtain Eq. ([3]). In doing so the expression of (fT3]) for 
p{z) has also been employed. 


Appendix D: Derivations for Eq. (1551) . and Eqs. (1621) 
and (16611 

For limiting large z, the leading-order contribution to 
the integral over k' on the rhs of Eq. ((54|) arises from the 
region of fc' —>■ fc, and hence we have 

4 K 

P,(z,z') = -^^ / dk \\cl>l{z)Mz') 

TT Zp(Z) Jo K 

4’k{z)(j}l{z'). (Dl) 

Similarly the contribution to the integral over k arises 
from the region of k ^ kp. Therefore 

p,(z,z') = / dkX^lMzW- 

TT Zp[Z) C Jo 

(D2) 
























13 


We then make the use of the expression (IT^ for p{z). 
Equation (1551) follows. 

We next present the derivation leading to Eqs. (l62ll 
and dnn). The emlicit expression for F{q) has been re¬ 
ported in Refs. [3, [1^ as 

F{q) = 7rXl0{X^ - A< - (?) 

-f6»(A> A< - q)e{q - X> + A<) 

x[7rA^ ^'{(A> - A^- (?} + Sx^{q) + Sx^{q)], 

(D3) 


Care therefore is needed in obtaining Eq. (lD9l) . We next 
list the following properties: 

A2>-X2^=A2<-X2^, (D12) 

and 

= ^5a<(<?) = -^Xl-Xl^. (D13) 

They are found to be useful in our further simplifying 
Eq. dml to 


where A> (A<) is the larger (smaller) one of A and A', 
respectively. In Eq. dm]), 


SM = - X.(A2 - Xlf/\ 


Sy{q) = Sx{q)\ 




and 


^ ^ ( A2 - A'\ 

^a(9) = :z{q + — -5 —); 
z q 


1 — A'^ 

Xy{q) = -iq-—^). 


(D4) 

(D5) 

(D6) 

(D7) 


The fact that the integral dqF{q) can be carried out 
analytically is the key point in obtaining the neat final 
expressions of (l62|) and (|M1) for bx{z, z'). To this end, we 
write it as 


/.oo rx/>‘l->-< 

/ dqF{q) = / dq-nX^ 

JO JO 

+ dq[SxM + SxA<i)]- (D8) 


/A>-A< 


The first integral on the rhs of the above equation is 
trivial. After performing a partial integration over q in 
the second one, we come to 

1-00 /'A>-l-A< o 

/ dqF{q) = - dq—[Sx^{q) + Sx^{q)]- 

Jo JA>-A< 

(D9) 

We caution the reader that the function Sx^{q) + Sx^{q) 
has a discontinuity at q = \/A> — A<: 


^a>(Ja2>-A2< ) + SxAJxi-x% ) = --x%, 


(DIO) 


^A> (J A2> - Al,) + 5., ( Ja| - A^ ) = -A^<. 


2 


(Dll) 


poo ^A>4-A< 

/ dqF{q) = 2 / dqqiXl - (D14) 

' 0 J A> —A<; 


Equation (ID14I) is rewritten trivially as 
ca/a?.-a 


dqq{X%-Xiy/^ 


/ dqF{q) = 2 

' 0 J A> —A< 

+ 2 f ^ dqq(Xl-XL Y^^. (D15) 

x/^> ~A< 


The rhs of Eq. (ID15I) can be further simplified via the 
following transform in the integrals: 


1 


A 


Under this transform, 


= —(A2<-A2 )i/2. 


(D16) 


=-A<vU^ 


q = y X% — A<a:2 — A< x/1 — x'^ 
in the first integral on the rhs of Eq. (ID15I) ; but 


(D17) 


^A< = A<x/l-a:2, 
q = y^A> — A<a;2 -|- A<\/l — 
in the second one. The resultant expression is 

'•1 A^+A|(1-2x2) 


(D18) 


pOC pi 

/ dqF{q) = 4A< / dxx 
Jo Jo 


^iXi-Xix^)il-x^) 

(D19) 


The preceding integral could be efficiently evaluated nu¬ 
merically. We here prefer to express it in terms of com¬ 
plete elliptic integrals. To this end, we cite the following 
relations (t > 1) [2^: 


dxx^ 


1 


x/{t^-X^){l-X^) 
= t 


= t 






(D20) 
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and 


1 


/ dxx^ —;= _ 

/o y'(t2_^2)(l_^2) 


1 

“ 3^ 
1 

^ 3^ 




(21^ + l)K ( - 2{t^ + 1)E 


, (D21) 


where F and E are the elliptic integrals of the first and 
the second kinds, respectively. Substituting these rela¬ 
tions into Eq. (ID19I) . one has 


noo 

/ dqF{q) = 
Jo 


4A> 


Ki^j(Al-Al) 


.(D22) 


Further substitution of Eq. (ID22[) into Eq. (EH) then 
leads to 


b^{z,z')= - 


A r^F nkp 


[^k,k'iz,z') + ^k',kiz^^')] 

A 


k(^)(A'2-A2) + e(0(A2+A'2) 

(D23) 

Equation (1621) follows from the symmetry property ^k,k' 
(z,z') = (z',z). 


We next make the variable transform y = A'^/A^ in 
Eq. (1^ to obtain, for z —>■ oo. 


pKp 

G{z,z') = (/'feF(z)|(^fep(z')P / dkX^(j)l{z) 

Jo 

f ^-c{kF-y/k%-X^v)z 

Jo ‘Jy'kj, - A2y 

x[K(v^)( 2 /-l) + E(V^)(y + l)]. (D24) 

Further transform x = {kp — k)z then yields 

G{z oo, z') = 


\<PkAz)V\^kAz)\\ (D25) 


where 


7 = 


x[K{Ay)iy-l) + EiAy){y + l)]. (D26) 


Finally we substitute Eq. (ID25I) into Eq. (l62)l and apply 
the identity: 

r dxx^Ge-<y+i)x ^ + 1)]-7/2 (D27) 

Jo ° 


to Eq. (jD26|) . The result is the expressions of (|55)) and 
(1571) . In doing so, we have also made the use of Eq. (1571) . 
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